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Abstract 

The decay width, forward-backward asymmetry and r lepton longitudinal and 
transversal polarization for the exclusive B — > Kt + t~ decay in a two Higgs doublet 
model are computed. It is shown that the forward-backward asymmetry and lon- 
gitudinal polarization of the r lepton are very effective tools for establishing new 
physics. 
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1 Introduction 



The analysis of flavor changing neutral current (FCNC) decays is one of the most promis- 
ing directions in particle physics, theoretical as well as experimental, as a potential testing 
ground for the Standard model (SM) and as regards to the physicists' endeavor to comply 
fully for establishing new physics beyond the SM Along these lines, the rare B meson 
decays which takes place via the FCNC, play an exceptional role. For example, an investi- 
gation of these rare decays opens the way for the possibility of a more precise determination 
of the Cabibbo-Kobayashi-Maskawa (CKM) matrix elements [Qj. 

Currently the main interest on the rare meson decays is focused on the decays for 
which the SM predicts the largest branching ratios that can be measurable in the near 
future. The rare B — > K£ + £~(£ = e, /i, r) processes are such decays. For these decays 
the experimental situation is quite promising with e + e~ and hadron colliders focusing 
only on the observation of exclusive modes with lepton pairs as the final states. The 
B — > K£ + £~ decay, which is described by b — > s£ + £~ transition at quark level, has been 
investigated extensively in both SM and two Higgs doublet model (2HDM) ||-[I6|. It 



is well known that in the 2HDM, the up type quarks acquire their masses from Yukawa 
couplings to the Higgs doublet H 2 (with the vacuum expectation value t> 2 ) and down type 
quarks and leptons acquire their masses from Yukawa couplings to the other Higgs doublet 
Hi (with the expectation value t>i). In 2HDM there exist five physical Higgs fields: neutral 
scalar H°, h°, neutral pseudoscalar A and charged Higgs bosons H ± . Such a model occurs 



as natural a feature of the supersymmetric models [T7|. In these models the interaction 



vertex of the Higgs boson and fermions depends on the ratio tan(3 = ^ which is a free 

parameter in the model. The constraints on tan/3 are usually obtained from B — B, K — K 
mixing, h — > 37 decay width, semileptonic decay b — > crv T and is given by [|TI| |19|j : 



(the lower bound m#+ > 200 GeV is obtained in ||19[). 

In all these studies the contributions from neutral Higgs boson exchange diagrams are 
neglected, since the lepton-lepton-Higgs vertices are proportional to the lepton mass. But 
for the b — > st + t~ decay the mass of the r lepton is not too small compared to the b 
quark mass, and hence one expects that the neutral Higgs boson exchange diagrams may 
contribute considerably to such channels. It has been pointed by Hewett [pj that the 
longitudinal polarization Pi of the final lepton is an important observable that may be 



accessible in the B — > Kr + T~dec&y mode. Recently it has been shown in [2TJ that the 
complementary information is contained in Pi, together with the two other orthogonal 
components of polarization (Pt is the component of of the polarization lying in the decay 
plane and Pn is the one that is normal to the decay plane). Both Pt and Pn are crucial 
for the t + t~ channel since they are proportional to j^. The b — > s£ + £~ transition 

contains three Wilson coefficients Cj, Cg^ and C10 in the SM. The different components of 
the polarization, i.e., Pi, Pt and P^, involve different combinations of Wilson coefficients 
CV, Cg^ and C10 (see below) and hence contain independent information. For this reason 
the polarization effects are thought to play an important role in further investigations of 
the structure of the SM and for establishing new physics beyond it. 
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The query for the calculation of the branching ratios and other observables requires 
the computation of the matrix element of the effective Hamiltonian responsible for the 
B — > Kt + t~ decay between B and K states. This problem is related to the non-perturbative 
sector of QCD and it can be solved only by means of a non-perturbative approach. 

These matrix elements have been investigated in the framework of different approaches 
such as chiral theory |22|, three point QCD sum rules method [p^j , relativistic quark model 
by the light-front formalism [Q, effective heavy quark theory |25| and light cone QCD 
sum rules |[26|| . The aim of the present work is to calculate these matrix elements in the 



light cone QCD sum rules in the framework of the 2HDM, taking into account the newly 
appearing operators, Cq^ and to study the forward-backward asymmetry and final lepton 
polarization for the exclusive B — > i^r + r~decay Taking into account the additional neutral 
Higgs boson exchange diagrams, the effective Hamiltonian is calculated in |27 



as 
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where the first set of operators in the curly brackets describe the effective Hamiltonian 
responsible for the b — > sl + l~ decay in the SM. Note that the contributions arising from 
the diagrams containing the charged Higgs bosons are taken into account by modifying the 
corresponding Wilson coefficients. These diagrams do not induce any additional operators. 
Their explicit forms and the corresponding Wilson coefficients Cj can be found in ||. The 
second set of operators in the brackets, whose explicit forms are presented in [p? 



come 



from the exchange of the neutral Higgs bosons. The corresponding Wilson coefficients are: 
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The QCD correction to the Wilson coefficients Ci(mw) and Cq^itiw) can be calculated 



using the renormalization group equations. In |27] it was shown that the operators Og and 
Oio do not mix with Qi (i = 1, . . . , 10), so that the Wilson coefficients Cg and Cio remain 
unchanged and their values are the same as in the SM. Their explicit forms can be found 
in |27|, where it is also shown that O7 can mix with Qi. But additional terms due to this 
mixing can safely be neglected since the corrections to the SM value of CV arising from 
these terms are less than 5% when tan/3 < 50. 

Moreover the operators Oj (i — 1, . . . , 10) and Qi (i = 3, . . . , 10) do not mix with Qi 
and Q2 and also there is no mixing between Qi and Qi- For this reason the evolutions 
of the coefficients Cq 1 and Cq 2 are controlled by the anomalous dimensions of Qi and Q2 
respectively: 

C Q Xm b ) = V^ Q/Po C Qi {m w ) , i = 1, 2, 
where 7q = —4 is the anomalous dimension of the operator S]J)r. 



Neglecting the strange quark mass, the matrix element for b — > sr + r decay is |27 



- 2CV^cv^(l + l5 )bfYr + C Ql s{\ + l5 )bfr + C Q2 s(l + 75)&r 75 r| , (8) 

where p 2 is the invariant dileptonic mass, the Wilson coefficients CV, Cg and C\q are ob- 
tained from their SM values by adding the contributions due to the charged Higgs bosons 
exchange diagrams. Note that this addition is performed at high scale, and then us- 
ing the renormalization group equations, the coefficients are calculated at lower m& scale. 
Coefficients Cq 1 and Cq 2 describe the neutral Higgs boson exchange diagrams' contribu- 
tions. Note that the coefficient Cg^(fi,p 2 ) = Cg(fi) + Y(fi, p 2 ), where the function Y 
contains the contributions from the one loop matrix element of the ||| ^8], p9|] . In addition 



to the short distance contributions, it is possible to take into account the long distance 
effects associated with real cc in the intermediate states, i.e., with the cascade process 
B — > KJ/i/)(i/j') — > K£ + £~. These contributions are taken into account by introducing a 
Breit-Wigner form of the resonance propagator and this procedure leads to an additional 
contribution to Cg^ of the form [ 1~0~[ , |30[ 
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From eq.(H) it is obvious that, in order to calculate the decay width and other observ- 
ables for the exclusive B — > K£ + £~ channel, the matrix elements (K |s7^(l — 7s)&| B), 
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(K \sia ia/ q 1/ (l + 7s)6| B), and (K \s(l + j 5 )b\ B) have to be calculated. These matrix el- 
ements can be parametrized in terms of the formfactors / + , /~ and fa in the following 
way: 

(K(p K )\s lfl (l- l5 )b\B(p B )) = ( P b+Pk)J + (p 2 )+pJ-(p 2 ) , (9) 
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m B + m K 

where p = p B — Pk is the momentum transfer. To be able to calculate the the matrix 
element (K \s(l + js)b\ B), we multiply both sides of eq.(|H]) by p M and use the equation of 
motion. Neglecting the mass of the strange quark, we get: 

1 
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Making use of eqs. 
i^T + r - decay: 

Ga 



1Q|) and ([Tl|) we obtain for the matrix element of the B 
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The formfactors f + {p 2 ), f~{p 2 ) and frip 2 ) are investigated in the light cone QCD sum 
rules framework and to a good accuracy their p 2 dependence are found to be representable 
in the following pole forms p6f : 

0.29 



which we will use in the numerical calculations. Using eq.(|T2"|) and performing summation 
over final lepton polarization, we get for the double differential decay rate: 
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where z = cos9 and 9 is the angle between the three momenta of the negatively charged 
lepton and the 5-meson in the CM frame of the final leptons, and v — Jl — is the lepton 



velocity Here A(l,r, s) is the usual triangle function and r = — s = As we have 
noted previously, the forward-backward asymmetry Ap B and the final lepton polarization 
involve different combinations of the Wilson coefficients Cj, Cg^ , Cio, Cq 1 , and Cq 2 and 
therefore each of them contains independent information. For this reason, here in what 
follows we study these quantities in more detail. 

The forward-backward asymmetry Ap B is defined as: 
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Note that in the SM, the forward-backward asymmetry is zero when the polarization 
of the final lepton is summed over. The reason is obvious: the hadronic current for the 
B — > K transition is a pure vector current. But forward backward asymmetry (or charge 
asymmetry) is non-zero only if there exists C-violating terms. In the 2HDM there is a C- 
violating term proportional to Fi (see eq.(|T2"|)), so that, A FB is nonzero and is proportional 
to the lepton mass. For small values of tanfl, the contributions from the neutral Higgs boson 
exchange diagrams are very small, and hence one expects that the value of A FB be small 
also. But for large tan/3, the contributions of the neural Higgs boson exchange diagrams 
become significant and we expect A FB to be large. The numerical analysis confirms these 
expectations (see numerical analysis section). 

Let us now discuss the lepton polarization effects. We define three orthogonal unit 
vectors: 
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where p\ and px are the three momenta of the i~ lepton and the K meson, respectively, 
in the center of mass of the £ + £~ system. The differential decay rate for any given spin 
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direction n of the I lepton, where n is a unit vector in the £ lepton rest frame, can be 
written as 



dT(n) 1 ( dT 
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(16) 



where the subscript "0" corresponds to the unpolarized case, and Pl, Pt, and P/v, which 
correspond to the longitudinal, transverse and normal components of the polarization vec- 
tor, respectively, are functions of p 2 . These components Pi (i = L, T, N) are defined 
as: 
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The calculations for the Pj's (i = L, T) lead to the following results: 
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The factor A in eqs. (|T^) and (|T^) can be obtained from eq.([15|) by an integration over z 
of the terms in the curly brackets. Note that the explicit form of the normal component 
Pjv °f the polarization vector of the l~ lepton is also calculated. However, an analysis of 
its behavior with respect to p 2 shows that numerically it is quite small, so that we do not 
present it. As a check of our results, when we equate Pi and P 2 to zero, i.e., neglect the 
contributions from the Higgs bosons, we obtain the results of 



2 Numerical Analysis 



The values of the main input parameters, which appear in the expression for the decay 
width are: m b = 5 GeV, m c = 1.4 GeV, m T = 1.78 GeV, = 0.105 GeV, Aq C d = 
225 MeV, = 5.28 GeV, and = 0.495 GeV. We use the pole form of the formfactors 
given in eq. (|I4]) . For B meson lifetime we take r{Bd) = 1.56 x 10~ 12 



SMi 
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to the leading logarithmic approximation 



The expression Cg^ for the b 



Cio = -4.642 . 
s transition in the next to leading order approximation is 



given as (see for example |32| ) 
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Cl ff (m b ) = 

Cl M (m b ) + C^(m b ) + 0.124tw(s) + g(m c , s) (3d + C 2 + 3C 3 + C 4 + 3C 5 + C 6 ) - 
- ±g(m q , s) (C 3 + 3C 4 ) - ~<7(m 6 , s) (4C 3 + 4C 4 + 3C 5 + C 6 ) + 

+ ~ ( 3C7 3 + ^4 + 3C 5 + C 6 ) , (20) 
y 

with 

Ci = -0.249 , C 2 = 1.108 , C 3 = 1.112 x 10~ 2 , C 4 = -2.569 x 10~ 2 , 
C 5 = 7.4 x 10" 3 , C 6 = -3.144 x 10" 2 , C$ M (m b ) = 4.227, 
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where m q = —j-, s = The explicit forms of (raw), Cg (m w ) and C^ (m w ) can 
be found in [KJ. 

In the above expression tu(s) represents the one gluon correction to the matrix element 
Og and its explicit form can be found in |12[], while the function g(m q ,s) arises from the 
one loop contributions of the four quark operators 0\ — O e (see for example [33[]), i.e., 
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+ 1 0(1 — y q ) (in — ^ ^ ~~ * 7r ) + ^(2/9 ~~ l)arctan 



- 1 . 



(21) 



where -u = and s 1 = 

In Table 1 we list the three different sets of values for the masses of the Higgs particles 
m h o, m H ±, m H o and m^o that we use throughout the numerical calculations. 







m H ± 






mass set-1 


80 GeV 


200 GeV 


150 GeV 


100 GeV 


mass set-2 


250 GeV 


300 GeV 


100 GeV 


350 GeV 


mass set-3 


100 GeV 


400 GeV 


200 GeV 


150 GeV 



Table 1: List of the values for the masses of the Higgs particles. 



7 



Further, we choose sina = and the numerical calculations for the branching ratio 
are all carried at three different values of tan/3, i.e., tan/3 = 1, tan/3 = 20 and tan/3 = 30. 

In Fig.l (a) we present the p 2 dependence of the differential branching ratio for the 
B d —> Kt + t~ decay with the long distance effects, for tan/3 = 1. The numbers 1, 2 and 3 
on each curve identify the mass set-1, 2 and 3 for the Higgs particles m#±, m h o, m H o and 
nrijfii respectively, as displayed in Table 1. The curve numbered as 4 is the one calculated in 
the Standard Model and is depicted for a comparison of the two models. Fig.l (b) and 
(c) are the similar graphs at tan/3 = 20 and tanj3 = 30, respectively. The sharp peaks in 
these figures are due to the long distance contributions. Note that the curves that represent 
the short distance contributions are not plotted in this set of figures. The reason for this 
is our observation of the fact that, in each respective mass set these curves overlap with 
and mimic the behavior of the curves representing the long distance effects at all points 
of p 2 , except at sharp peaks. It is also observed that the spectrum of the invariant mass 
distribution is slightly asymmetric. 

In Fig. 2 (a) and (b) we plot the dependence of the forward-backward asymmetry Afb 
on p 2 for the Bd — > Kt + t~ decay, with and without the long-distance effects at different 
values of tan(3. The lines numbered as 1, 3 and 5 represent the long distance effects for the 
mass sets 1, 2 and 3, respectively, while the lines with numbers 2, 4 and 6 are depicted for 
the short distance effects. From these figures we see that Afb is negative for all values of 
p 2 except in the tp' resonance region and it is sensitive to the value of tan(3. For tan/3 = 1 
Afb is quite small, so that we do not present it here. 

In Fig. 3 (a), (b) and (c) we present the p 2 dependence of the longitudinal polarization 
of the final lepton P L without the long distance effects at tan[3 = 1, tan[3 = 30 and 
tan/3 = 50, respectively. The lines numbered as 1, 2 and 3 represent the mass set-1, 2 and 
3, respectively. The line with number 4 is the one presented in |26| for the Standard Model 
calculations. As is obvious from these figures if we exclude the resonance mass region of 
ip' , Pl is negative for all values of p 2 . Note that the behavior of the Pl curves are sensitive 
to the different choices of the value of tan/3. 

In Fig. 4 (a), (b) and (c) we present the p 2 dependence of the transversal polarization P? 
of the r lepton which lies in the decay plane, without the long distance effects, at tan(3 = 1, 
tan(3 = 30 and at tan/3 = 50, respectively. The lines numbered as 1, 2 and 3 represent the 
mass set-1, 2 and 3, in the respective order. From these figures it follows that at tan/3 = 1 
Pt is positive at all values of p 2 . On the other hand, for tan(3 = 30 and tan(3 = 50, Pt is 
positive near the threshold region while it becomes negative far from the threshold region. 
Therefore the determination of the sign of Pt in the future experiments is a very important 
issue and can provide a direct information for the establishment of new physics. 

For completeness in Figs. 5-7 we present the results for the B — > Kfi + u~dec&y. It is clear 
that in this case the neutral Higgs exchange diagram contributions are quite small and the 
deviation from the SM prediction is due to the charged Higgs boson exchange diagrams. 

In Fig. 5 we present the differential branching ratio versus p 2 for B — > K{i + u~with and 
without the long distance effects. The sharp peaks in these figures are due to the long 
distance contributions, as is the case for Fig.l. Fig. 5 (a) is for the mass set-1, (b) is for 
the mass set-2 and (c) is presented for the mass set-3. The lines numbered as 1, 3 and 
5 represent the long distance effects for the mass sets 1, 2 and 3, respectively, while the 
numbers 2, 4 and 6 are presented for the short distance effects. The abbreviation SM on 
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the bottom pair of curves (5 and 6) stands for the Standard Model results as depicted in 
2~6fl . The branching ratio for the B — > fT/i + /i~decay is not as sensitive to the value of tanfi 



as for the B — > i^r + r _ decay. 

The behavior of the longitudinal polarization Pl with changing p 2 for the B —>■ Kfi + u~ 
decay is presented in Fig. 6, with and without the long distance effects. Curves 1, 2 are 
evaluated numerically at tan(3 = 1 while those 3 and 4 represents the behavior of the 
longitudinal polarization calculated at tan/3 = 50. The last two curves 5 and 6 depict the 
Standard Model results f2"6fl . The ordering of the figures and of the numbering of lines for 
the short and long distance effects are exactly the same as explained for Fig.5. Without 
the long distance effects Pl is always negative. 

In Fig. 7 (a), (b) and (c) we present the transversal polarization Pt for the B — > 
K n + \i~ decay as a function of p 2 , for mass set-1, mass set-2 and mass set-3, respectively. 
While at tan[3 = 1 Pt is always positive, at tan(3 = 30 and tan[3 = 50 it changes sign. 
Thus the investigation of the sign of Pt can be an effective tool in search for new physics. 

Since, as we have already noted, Apg, Pl, and Pt contain independent information, 
their investigation in the future experiments will be quite an efficient tool for establishing 
new physics. 

In Table 2 we present the values of the branching ratios for the Bd — > Kt + t~ decay. 
After integrating over p 2 we get for the branching ratios for the Bd — > Kt + t~ decay, with 
the long distance contributions: 





tanf3 = 1 


tanp = 20 


tanf3 = 30 


mass set-1 


3.47 x 10~ 7 


2.68 x 10~ 7 


4.12 x 10~ 7 


mass set-2 


3.15 x 10~ 7 


2.64 x 10~ 7 


3.71 x 10~ 7 


mass set-3 


2.96 x 10~ 7 


3.01 x 10~ 7 


5.43 x 10~ 7 



Table 2: Branching ratios for the exclusive Bd — > Kt + t decay, with the long distance 
contributions. 



The ratio of the exclusive and inclusive channels is defined as: 

B(B d -> Kt + t-) 
B{b -> st + t-) ' 
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In the SM this ratio is given as R = 0.07±0.02 when B(b st + t~) = (2.6±0.5) x 10~ 7 
33fl . In Table 3 we display the results that we have calculated for our case, using the values 



given in Table 2. 





tanf3 = 1 


tanf3 = 20 


tanft = 30 


mass set-1 


1.33 


1.03 


1.58 


mass set-2 


1.21 


1.02 


1.43 


mass set-3 


1.14 


1.16 


2.09 



Table 3: The ratio R of the exclusive B(B^ — > Kt + t ) and SM value of the inclusive 
B{b — > st + t~) channels. 



where we have used the SM value for the inclusive B(b — > st + t~). 

For completeness we present in Table 4 the branching ratio with respect to the inclusive 
B(b — > st + t~) decay in 2HDM with the long distance contributions. The expression for 
the inclusive B(b — > st + t~) decay is given in [27 . 





tanf3 = 1 


tanp = 20 


tanf3 = 30 


mass set-1 


18.90 x 10~ 7 


17.97 x 10~ 7 


18.44 x 10~ 7 


mass set-2 


18.61 x 10~ 7 


17.98 x 10~ 7 


18.34 x 10~ 7 


mass set-3 


18.44 x 10~ 7 


18.13 x 10" 7 


18.98 x 10~ 7 



Table 4: Branching ratios for the inclusive Bd — > Kt + t decay in the 2HDM, with the 
long distance contributions. 
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Finally in Table 5, we display the results for the ratio R of the inclusive B(b — > st + t ) 
decay calculated in the 2HDM and SM value of the inclusive B(b — > st + t~) channel. 





tanj3 = 1 


tan(3 = 20 


tanf3 = 30 


mass set-1 


7.27 


6.91 


7.09 


mass set-2 


7.16 


6.92 


7.05 


mass set-3 


7.09 


6.97 


7.30 



Table 5: The ratio R of the inclusive B(b -> st + t~) decay in 2HDM and SM. 



In conclusion, we calculate the rare B — > K£ + £~ decay in 2HDM. It is observed that the 
forward-backward asymmetry A FB , the longitudinal polarization P L and the transversal 
polarization P T of the charged final lepton are very sensitive to the variations in tanf3. 
Therefore, in search of new physics their experimental investigation can serve as the crucial 
test. 
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Figure Captions 



1. Invariant mass distribution for the B — > Kt + t~ decay with the long distance effects. 
In Fig.l (a), (b) and (c) the numbers 1, 2 and 3 on each curve identify the mass set-1, 2 
and 3, respectively, for the Higgs particles. The curve numbered as 4 is the one calculated 
in the Standard Model. The sharp peaks in these figures are all due to the long distance 
contributions. 

2. The dependence of the forward-backward asymmetry Afb on p 2 for the decay B — > 
Kt + t~ . The lines numbered as 1,3 and 5 represent the long distance effects for the mass 
sets 1, 2 and 3, respectively, while the lines with numbers 2, 4 and 6 are chosen for the 
short distance effects. 

3. The dependence of the longitudinal polarization, Pl, on p 2 for the B — > Kt + t~ without 
the long distance effects. The lines numbered as 1, 2 and 3 represent the mass set-1, 2 
and 3, respectively. The line with number 4 is the one presented for the Standard Model 
calculations. 

4. Transversal polarization asymmetry, Pt, for the B — > 7^T + T~decay as a function of 
p 2 . In this set of figures the lines numbered as 1, 2 and 3 represent the mass set-1, 2 and 
3, respectively. 

5. The same as in Fig.l, but for the B — > Kfi + fi~ decay. Fig. 5 (a) is for mass set-1, 
(b) is for mass set-2 and (c) is presented for mass set-3. The lines numbered as 1, 3 and 
5 represent the long distance effects for the mass sets 1, 2 and 3, respectively, while the 
numbers 2, 4 and 6 are presented for the short distance effects. The abbreviation SM on 
the bottom pair of curves (5 and 6) stands for the Standard Model results. 

6. The same as in Fig. 3 but for the B — > Kfj, + fi~ decay. The ordering of the figures 
and of the numbering of lines are exactly the same as explained for Fig. 5. 

7. The same as in Fig. 4, but for the B — > Kfj, + fi~ decay. Figures (a), (b) and (c) are 
presented for the mass set-1, mass set-2 and mass set-3, respectively. 
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